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Abstract 

One dimensional systems are under intense investigation, both from theoretical and experimental points of view, since 
they have rather peculiar characteristics which are of both conceptual and technological interest. We analyze the depen- 
dence of the behaviour of one dimensional, time reversal invariant, nonequilibrium systems on the parameters defining 
their microscopic dynamics. In particular, we consider chains of identical oscillators interacting via hard core elastic col- 
lisions and harmonic potentials, driven by boundary Nose-Hoover thermostats. Their behaviour mirrors qualitatively that 
of stochastically driven systems, showing that anomalous properties are typical of physics in one dimension. Chaos, by 
itslef, does not lead to standard behaviour, since it does not guarantee local thermodynamic equilibrium. A linear relation 
is found between density fluctuations and temperature profiles. This link and the temporal asymmetry of fluctuations of 
the main observables are robust against modifications of thermostat parameters and against perturbations of the dynamics. 

1 Introduction 

Understanding the fluctuation properties of nonequilibrium phenomena is one of the major tasks of modern statistical 
physics Ill|2][3][4][5]|6]|7]E]|9]. T° tn i s purpose, various generalizations of Onsager-Machlup theory [ 10] have been pro- 
posed, like those of Refs.|(8] |9) . These papers have also motivated various works on deterministic nonequilibrium particle 
systems IfTTI fT2l [T3l fT4l IT31 meant to investigate the relationship between stochastic and deterministic models of nonequi- 
librium physics. One basic tenet of classical statistical mechanics, indeed, maintains that the stochastic description is but 
a reduced (mesoscopic) representation of the deterministic (microscopic) description of a given system of interest. Nev- 
ertheless, the relation between the two representations is far from fully understood. For instance, nonequilibrium steady 
states require thermostats, and one expects the state of a macroscopic system to be unaffected by the details of the ther- 
mostatting mechanism. Thus, even in mathematical models of nonequilibrium steady states, one would like stochastic and 
deterministic thermostats to lead to practically equivalent representations of the same phenomenon^ However, different 
theoretical models may capture different aspects of a physical process, and a complete equivalence should not be ex- 
pected, except in some limiting situation, as postulated by various equivalence principles lfT7l[T^[T^l^l2Tll22ll23l . The 
presence of local thermodynamic equilibrium fosters the equivalence. Indeed, the thermodynamic relations are so largely 
independent of the details of the microscopic dynamics, hence so generally valid, precisely because of local equilibrium, 
cf. Section 4 below, Refs.[2 24 25 26]|27][28] and references therein. In case local equilibrium cannot be established, 
one may still wish to organize the different cases in homogeneous groups, or universality classes 11291 l30l . to help our 
understanding of nonequilibrium phenomena. 

1 Rieder, Lebowitz and Lieb, speaking of nonequilibrium systems due to stochastic boundary reservoirs at different temperatures, express this idea 
as follows: "the properties of a 'long' metal bar should not depend on whether its ends are in contact with water or with wine 'heat reservoirs' 
at temperature T\ and T2" 1161 . In other words, the temperature should be the only thermodynamically relevant property of a thermostat, in the 
nonequilibrium steady state. 
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It is, then, interesting to point out the conditions under which the asymmetries of fluctuations characterizing the (in- 
trinsically irreversible) stochastic evolutions |[8][3T) are present in the dynamics of time reversal invariant particle systems. 
These asymmetries may indeed be experimentally and numerically observed, and are thought to be responsible for the 
irreversibility of macroscopic phenomena BTI . In Ref. IfTTI . the fluctuations of the current of the nonequilibrium Lorentz 
gas, with large numbers N of noninteracting particles, were found to be time-symmetric, despite the clearly irreversible 
behaviour of the system. As explained in Ref. [14], this was due to the lack of interactions among the particles, hence 
to the lack of correlations among them [15), revealed by the large system limit@ On the other hand, the absence of 
interactions prevents the onset of local thermodynamic equilibrium, despite the convergence to a given steady state. In 
Refs. |[T2l[T3l[T4l[l"5l , it is therefore argued that systems of properly interacting particles should typically have asymmetric 
fluctuation paths, as predicted for stochastic systems. More precisely, Ref.[ 14 1 states that temporally symmetric fluctua- 
tion paths of nonequilibrium time reversal invariant deterministic dynamics require the untypical condition of vanishing 
correlations; generically, fluctuation paths should then be asymmetric. Hence, deterministic reversible and stochastic 
interacting particle systems are equivalent, i.e. belong to the same universality class, from this point of view. 

In the present paper, we continue the investigation of these issues, by looking at variations of the nonequilibrium 
Fermi-Pasta-Ulam (FPU) model introduced by Lepri, Livi and Politi [ 33 1 . This model consists of N anharmonic oscillators 
of same mass m and positions xj, j = 1, ...,N, interacting via the nearest neighbours FPU-yS potential [34|: 



V(q) = k 2 ^+/3^, (1) 

while the oscillators with /' = 1 and j = N further interact with deterministic thermostats known as Nose-Hoover "thermal 
baths", and with still walls. The internal energy of the system is then given by: 

N 1 

H=Y,2 m ^ +v(q j +i - q ' y (2) 



The equations of motion take the form 



mqi = F(q x - q ) - F(q 2 - q x ) - &q\ , 

mqj = F(qj - q M ) - F(q j+l - qj) , for j = 2 N- 1 (3) 

mq N = F(q N - q N _ x ) - F(q N+l - q N ) - % r q N , 

where qo = qu+i = represent the walls and m is the mass of one particle, which we take to be unitary, so that velocities 
qj and momenta p, represent the same quantity. Moreover, the interparticle forces are obtained differentiating the poten- 
tial, F(q) — -V'(q), and the Nose-Hoover thermostats, at "temperatures" T r and 77 with response times 6 r and 8(, are 
implemented by the variables & , £ r , which obey 



. _ 1 



1 



-1-1 

Te 



(4) 



The nonequilibrium FPU-j6 model is time reversal invariant, but dissipative, in the sense that the time average of the 
divergence of the equations of motion, -((£ r ) + ((e)), is negative. Therefore, the sum of the Lyapunov exponents is 
negative, the phase space volumes contract and the system approaches in time a nonequilibrium steady state, characterized 
by a singular phase space probability distribution. 

In the theoretical calculations, the motion of particles is commonly referred to the quantity a, which is the mechanical 
equilibrium distance between two nearest neighbours [25 1. Then, the mean displacement of the position of particle j from 
its equilibrium position ja, is often assumed to be small, in order to define microscopically the physical observables. It 
will be pointed out that this is not appropriate, except for particularly small driving forces. 



"Asymmetric fluctuations are found at small N, since the averages computed over small numbers of particles do not accurately reproduce the 
correlation functions [ 32 ]. 
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A quantity of interest is the "local virial" ll33l . For particles in the bulk (z = 2, ...,N - 1), this is the time average 
of the product of the displacement of each particle, times the net force acting on it, (<7,(^/+i _ Ff)), where we have set 
Fj - F(qi - qi-\) for simplicity. If the anharmonic part of V is replaced by hard core elastic collisions, the interaction 
between particles is purely harmonic, and one has Fj - F,-+i = fc 2 (g,-i - 2qi + qi + \). Therefore, the local virial of particles 
in the bulk is expressed by 

k 2 (-qi (qi-i - 2qt + <7,+i)> , i = 2, ...,N - 1 

while for i = \,N it reads 

(qi - k 2 (2qi - q 2 ))) , (q N (~^RqN ~ k 2 (2q N - q N -i))} ■ 

An important feature of this model is that its long wavelength Fourier modes, representing the slow relaxational dynamics, 
may be considered as approximately independent of the short wavelength modes, representing the fast dynamics [25], 
i.e. representing some sort of noise added to the slow dynamics^] In the slow dynamics, one may further separate a 
conservative non vanishing harmonic part, which cannot contribute to transport phenomena, from a mode interaction part 
[25, 36 , 37 1 . From this standpoint, the nonequilibrium FPU system is similar to the nonequilibrium Lorentz gas, as argued 
in [[14 1, but the presence of nonvanishing correlations, even if decaying in time, provides a mechanism for asymmetric 
fluctuation paths, within the corresponding decorrelation time scales, Refs. |fT3l [T4 lFI 

In this paper, we consider different values of k 2 and a different kind of anharmonicity: that provided by hard core 
elastic collisions. We also consider the limiting cases in which either the harmonic or the anharmonic interactions vanish. 
These limiting cases share indeed some peculiarities of noninteracting particle systems, and constitute good candidates 
for the study of the onset of temporally symmetric fluctuations |[T4l [151 . The fact that all particles are identical, further 
means that there is no disorder in our chains, hence no reasons for chaos, mixing or ergodicity, except those that may be 
contributed by the boundaries. We will see that nonequlibrium boundary conditions make correlations develop even in the 
bulk of such particles chains, and lead to asymmetric fluctuations. The results presented here lead to the following main 
conclusions: 

a) The kinetic temperature profiles are linearly related to the deviations of the mean positions of the particles from 
their equilibrium values, hence to the deviations from uniformity of the density profiles (Section 2): 

Ti=Mxi+i-Xi)+p2. (5) 

This relation and the temporal asymmetry of fluctuations are robust against modifications of the parameters of the 
thermostats and of the interaction potentials, besides being robust against stochastic perturbations of the determin- 
istic dynamics (Appendix B). 

b) Comparison between our results for deterministic, time reversal invariant, dynamics and those found in the literature 
for chains with stochastic thermostats, shows that the behaviours concerning the two kinds of thermostats are only 
qualitatively similar. The two kinds of thermostats do not belong to the same universality class, except in a coarse 
sense; equivalence is then only expected in the equilibrium limit. 

Other relevant conclusions include: 

c) The deviations of the mean positions of the particles from their equilibrium values are large, except in the purely 
harmonic case which, however, does not sustain a temperature gradient. Hence the microscopic definitions of 
nonequilibrium thermodynamic observables, based on small deviations, are generally impaired (Section 2). 

d) As in stochastically thermostatted systems, investigated by other authors, the properties of steady states depend 
substantially on the microscopic details of the dynamics, showing that the absence of genuine local thermodynamic 
equilibrium is typical of physics in 1-d (Section 2). 

3 See, e.g. Ref.| 35 1 for a discussion of the relation between fast and slow variables. 

4 The role of correlations for temporal asymmetries has been emphasized in Ref. [151 , although in the different framework of homogeneously driven 
and thermostatted systems. 
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e) As in the stochastic case, temporal asymmetries are ubiquitous in deterministic, time reversal invariant, nonequilib- 
rium particle systems, as long as correlations among particles are present. Such asymmetries are not restricted to 
large deviations in the large system limit and do not require local thermodynamic equilibrium (Section 3). 

f) At variance with common expectations, Nose-Hoover thermostats alter the bulk behaviour, even when the bulk 
resembles a noninteracting particle system, for both small and large temperature gradients. This is a manifestation 
of the nonlocality generically observed in nonequilibrium steady states |8, 38, [39j |40) (Section 4). 

g) Chaos by itself is not sufficient to establish standard behaviour in 1-d systems, besides not even being necessary 
113311411 1421 1431 (Section 5). 

In Section 2, we consider hard spheres on a line, connected in pairs by harmonic forces. In Section 3, the temporal 
symmetries of the fluctuations of this model are analyzed. In Section 4, the limiting cases of purely harmonic forces and 
purely elastic collisions are considered. In section 5, the role of chaos is investigated. Section 6 is devoted to concluding 
remarks. Appendix A quantifies the deviations of the temperature profiles from a theoretical curve. Appendix B gives 
results about stochatstically perturbed dynamics. 
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2 Hard core and harmonic potentials 

In the present case, the harmonic part of the potential of Eq.(fl]i is retained, while the term proportional to (3 is replaced by 
hard core elastic collisions between particles of radius r — 1 . This model shares some features of non-interacting particles 
systems. Indeed, isolated one-dimensional systems of elastic particles of equal mass and size exchange their momenta 
in such a way that the overall motion is equivalent to that of a system of non-interacting particles, in which each particle 
preserves its momentum [44 1. Moreover, in isolated systems, the harmonic potential should not constitute an important 
source of correlations, since the oscillations it entails amount to independent normal modes. Nevertheless, two features 
distinguish our model from systems of non-interacting particles: 

• the elongation of the harmonic springs discriminates the case in which particles bounce back at collisions, from the 
case in which they pass through one another; 

• the presence of the thermostats at the boundaries of the chain spoils the normal modes decomposition of the purely 
harmonic motion. 

The question, which we have investigated numerically, is how crucially these facts affect the behaviour of the system. To 
verify that the steady state is indeed achieved, we have computed the local "temperatures", (pif), and the virial profiles, 
at different times along a simulation of length ? max • We have found that the temperature profile rapidly approaches its 
asymptotic form (a quarter of our ? max suffices for good accuracy), while the virial expression takes much longer to settle 
down, cf . Figure 1 for the case with N = 400, k 2 = 1 , 8? = 8 r — 1 . 

Similar profiles have been found in numerous one-dimensional nonequilibrium chains of oscillators, with anhar- 
monic interaction potentials, like the FPU-yS chains, [25 45 46 1. In particular, assuming that dynamical chaos practically 
amounts to some degree of stochasticity in the evolution, our model should behave similarly to that of Ref. 11271 . which 
has stochastic boundary thermostats, harmonic interactions and stochastic exchanges of momenta between nearest neigh- 
bours, meant to mimic (momentum and energy preserving) hard core interactions. In a suitable continuum limit, the 
corresponding energy profile takes the form [ 27 1 : 



where £ is the Riemann ^"-function, and the space variable is rescaled so that the left end of the chain corresponds to 
x = — 1 while the right end corresponds to x = 1. 

However, as shown by the left panel of Figure 2, and deduced from a comparison of our Figure 1 with Figure 2 of l27l . 
the temperature profiles obtained by Lepri, Mejfa-Monasterio and Politi are not completely equivalent to ours. In the first 
place, the slope of profile © has square root divergences at the boundaries, which are less steep than ours. Secondly, the 
profiles of Ref. [27 1 are odd with respect to the centre of the chain, even at finite N, while our profiles are not. One further 
difference is the simple k 2 dependence of the finite-A^ temperature profiles of 11271 . as opposed to the strongly irregular 
dependence of our profiles, cf. central panel of Figure 2. To investigate this question, we have analyzed the case with 
k 2 = 1, T { = 320, T r = 20,6 { = 6 r = l, for N = 100, 200, 250, 300, 400. Differently from |27), we found no evidence of 
the Af~'/ 3 rate of convergence towards the asymptotic profile T(x) and, in most of the chain, the difference between our 
profiles and T(x) was observed to grow with N. Therefore, as shown in Appendix A, either a scaling regime towards T(x) 
sets in at N larger than 400, or our asymptotic profiles will not correspond to those of 11271 . On the other hand, Appendix 
A shows that our profiles converge rapidly to a given asymptotic profile, which is not odd. This may be related to the fact 
that, differently from the stochastic case, the energy exchange between Nose-Hoover boundary thermostats and 1-d (or 
almost 1-d) systems of particles is more efficient at the hot side than at the cold side, cf. Ref. flTl and Section|4]below0 

5 In Ref. 1471 it was observed that small temperature gradients allow the cold thermostat to behave as "thermodynamically" as the hot thermostat. 
The boundary thermostats of 1-d chains act only on momentum degrees of freedom, making it difficult, especially at large temperature gradients, for 
local equipartition to be established, except in the presence of mass diffusion |48 1. However, the difference between T and the profiles obtained by 
Nose-Hoover boundary thermostats persists at small gradients, as further evidenced in Sectionf4] 




(6) 
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Figure 1. (Color online) Kinetic temperature (left panel) and local virial (right panel) profiles for different fractions (1/4, 1/2, 3/4, 
1) of the simulation length /max = 2.5 • 10 s time steps of size h = 1CT 4 , with N = 400, k 2 = 1, ( = r = 1. The (p 2 ) profile reaches 
rather rapidly its asymptotic form (the curves for the four simulation lengths are indistinguishable), while the virial expression takes 
longer to converge. Because of hard core interactions, kinetic temperature and virial do not coincide 1491 . 



Despite these facts, the temperature profiles concerning nonequilibrium chains of oscillators typically enjoy the same 
qualitative behaviour, consisting of steep curves at the boundaries, due to contact resistance, interpolated by almost straight 
lines in the bulk. The nature of this behaviour is far from obvious, cf. the discussion in Ref. lt27l |46l . In any event, the 
strong and irregular dependence of the steady states of both deterministic and stochastic 1-d systems, on microscopic 
details of the dynamics, as well as the violations of Fourier Law, reveal the absence of genuine local thermodynamic 
equilibrium and of diffusion, see e.g. Refs. lfT6l l25, 26, 28 45]. 

For certain observables, this dependence on the microscopic mechanisms, hence the lack of local thermodynamic 
equilibrium, persists in the large /V limit, lfl6l l27l 1281 . Similarly, completely different situations are determined by the 
boundary conditions, such as the divergence of the conductivity of disordered harmonic chains with free boundaries, in 
the thermodynamic limit, as opposed to its vanishing trend in chains with fixed boundaries [28 1 30]. Moreover, the scaling 
behavior of the conductivity with system size crucially depends on the spectral properties of the heat baths l50ll . 
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Figure 2. (Color online) Comparison between the temperature profile T(x), eq. {6]l, of Ref.| 27 1 (continuous line) and those of our 
deterministic models, for N = 100, TV = 320, T r = 20 and decreasing k 2 (left panel). Dependence of the deterministic profile on k 2 
(central panel). For k 2 — » 0, it tends to the k 2 = 0-profile in an irregular fashion. Arrows indicate the values of k 2 . Deterministic 
profiles as functions of T ( , with k 2 = 0. 1 and T r = 20 (right panel). 
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Rieder, Lebowitz and Lieb interpreted this kind of results as "unphysical" 1161 . Indeed, were 1-d chains of oscillators 
to represent macroscopic objects enjoying thermodynamic properties, the thermostat dependence of the heat fluxes of 
|[T()1 , or of our temperature profiles, would be unrealistic. However, recent theoretical and technological developments 
allow different interpretations. The study of nearly 1-d systems, especially with non-macroscopic numbers of elemen- 
tary constituents, indicates that standard thermodynamic relations typically fail to describe their behaviour, although a 
complete understanding of their properties is currently missing, see e.g. QUI l46l ISTl l52l and references therein. Thus, 
a possible absence of local thermodynamic equilibrium does not need to be unphysical. Indeed, local thermodynamic 
equilibrium requires a sufficiently fast decay of correlations of all relevant observables, which would lead to diffusive 
transport. But this is seldom afforded by (quasi-) 1-d systems ||26l l53l FI which include modern technological artifacts 
as well as structures found in nature. It is also known that 1-d systems are affected by peculiarities such as the cumu- 
lative O(N) fluctuations, about the particles equilibrium positions |54|. These frustrate the microscopic definitions of 
observables based on the assumption of small fluctuations and, more importantly, are at odds with the properties of solids. 
Furthermore, we have found that temperature gradients induce large displacements of the average equilibrium positions, 
making the distribution of matter inhomogeneus. For three different values of k 2 , the left panel of Figure 3 portrays these 
average deviations, {qi) - (xj - ai), for i = 1, ...,N = 100. The right panel of the figure portrays the quantity (xt+\ - x,), 
which is related to the inverse of the density of particles. One observes that nonequilibrium boundary conditions lead to a 
shift of all particles towards the cold side of the chain, hence to a density gradient. 

Denoting by (•)# the averages computed with elastic constant equal to k 2 , closer examination shows that: 

• all {qi)p are positive and, in the bulk, are approximated by a parabola with maximum at i > N/2 for k 2 + 0, while 
the maximum of (g,)o occurs at i = N/2; 

• consistently, {qN)k 2 > (li)k 2 f° r k 2 + 0, while (^)o = (qi)o', 

• (1n)ii 2 = {<1n)p f° r a H pairs k 2 , k 2 , while (qi)^- = (qi)p only if both k 2 , k 2 + 0; 

• the dependence of (q,)^ on k 2 is not monotonic, like the temperature profile is not. 

Therefore, the systems with pure hard sphere interactions, whose bulk corresponds to noninteracting particles, behave 
differently from the systems with both hard core and harmonic interactions. Indeed, while the temperature profiles at 
k 2 — do not substantially differ from those at k 2 + 0, the net energy transfer with k 2 — is singular, since it is equivalent 
to that of a single particle, bouncing back and forth between the hot and cold walls, independently of the value of N. 




Figure 3. (Color online) Mean displacement {qt) k 2 from the equilibrium centres of oscillation of the particles, ai (left panel). Mean 
coordinates difference {x i+ \ — x ; ), for several vales of k 2 (right panel). The higher this difference, the lower the particles density. In both 
panels N = 100, T t = 320, T r = 20. 

6 For instance, certain correlations never decay in the narrow channels of Ref. 1531 and in our chains of hard core particles, since the particles order is 
preserved in time. While this is not a problem per se -order is preserved in 3-d crystals as well- it violates the hypothesis of molecular chaos and bears 
substantial consequences on transport phenomena. 
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Recalling that a — 5, in our calculations, these plots show that the steady state (average) particle distribution is not 
uniform and that the deviations of particles positions from their equilibrium values are large and correlated to the kinetic 
temperature profiles. Indeed, (jCj+i - which, apart from an unessential additive constant, is the discretized derivative of 
(qi), qualitatively approximates the typical temperature field (compare e.g. the right panel of Figure 3 with Figures 1 and 
2), hence a linear relation between the two quantities may be surmised: 

T, = filial -*)+ Pi- (7) 
The good agreement between T, and the kinetic temperature profile (p?) is demonstrated by Fig. 4. 
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Figure 4. (Color online) Comparison between temperature profiles (red circles) and rescaled average displacements of neighbouring 
particles Pi{x i+ \ - xi) + Pi (blue stars), for N = 100, TV = 320, T, = 20. The constants /?i,/?2 are obtained by a least-square fit of the 
data of Fig. 3. The values of the elastic constant are: k 2 = (left panel), k 2 = 0.1 (central panel), k 2 = 1 (right panel). 

This relation between temperature profiles and rescaled average displacements of neighbouring particles is robust against 
modifications of all the parameters of the dynamics, including thermostat parameters and interaction potentials, besides 
being robust against (small and large) stochastic perturbations of the dynamics, cf. Appendix B. Apart from the temporal 
asymmetries of the fluctuations of the main observables, considered in the next section, this is the only result which does 
not show a delicate dependence on the details of the microscopic dynamics. 
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3 Temporal symmetry of fluctuation paths 



In the literature, various notions of fluctuation path have been investigated. Given the substantial equivalence of the results 
based on such different notions, we adopt the first definition of fluctuation-relaxation path of Ref. [HI, denoted by FR1 Q 
Denote by ^# the phase space of our system, by S' : ^ — > the time evolution operator, and by X : ^# — > R 
an observable of interest. Consider an initial phase T e ./#, in the support of the steady state phase space probability 
distribution, and denote by X, the quantity X(S T), under the assumption that all, but a set of vanishing probability, such 
initial conditions enjoy the same statistics. Choose a fluctuation value ,^{X). The FR1 fluctuation path is defined as 
follows: 

Definition. Assume that X; — 3?(X). Then, for any ?o, t > 0, the FR1 fluctuation path of duration 2to based at i is the 
curve 

\X i+T : re [-f ,f ]} 

in the (r, X)-plane. 

As in Ref. [ 14], the symmetry properties of the fluctuation paths of the observable X are here assessed setting a threshold 
,^{X) - E(X) + 3cr(X), three standard deviations above the mean E(X), while the width of the observation time interval 
is taken to be 2fo = 2. Then, all such time intervals are translated by -t, so that the centers of all fluctuations coincide at 
t = 0. 

The first fluctuating observable we consider is the local dimensionless density of particles in the center of the chain, 
defined in a box of size L, ideally with 1 « I « W: 

p-po , x/ = minx;, x F =maxxj, (8) 

Xf — Xj jeL jeL 

where po = L/(L— l)a is the equilibrium density. 

For a family of n (time translated) fluctuations, {p w (f)> t £ [— fo, fo]} _ t > we obtained the average path p T subdividing 
[-to, to] in bo = 100 bins, and by computing 

-, n t 

p T = rr y y p s (t)%(t) (9) 
BJB ( T )^2-ji 

where r labels a bin, and %-(f) = 1 if t belongs to the r-th bin, while 1£ T (t) = if it does not. The quantity m(j) is the 
number of integration time steps which make a bin@ 

Figure 5 represents the normalized average paths, R T = (p T - minp T )/(maxp T - minp T ), which start at zero in the first 
bin and reach 1, their maximum value, in the central bin. The result is that the average path is temporally asymmetric, with 
growth steeper than relaxation, and that, differently from the observables considered in Section 2, it is not appreciably 
affected by k 2 . 

7 See fT3. 14] for a discussion of the ambiguities in the definitions of fluctuation paths in deterministic dynamics. Despite such ambiguities, the dif- 
ferent notions of fluctuation path have led to analogous conclusions and can be used interchangeably to asses the symmetry properties of nonequilibrium 
fluctuations. Furthermore, the different notions coincide in the large N limit, cf. 1 12. 13 14 15]. 

8 In our simulations, the integration time step equals 10~ 3 for small N, and 10~ 4 for large N, so that the numerical errors are of the same order. 
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Figure 5. (Color online) Average fluctuations of density in the center of the chain, for different chain lengths, N = 100 (left panel), 
N = 200 (central panel), N = 400 (right panel) and different values of the elastic constant, k 2 . The fluctuation paths are temporally 
asymmetric and not sensibly affected by the value of k 2 . 

This is consistent with Refs. 01411151 . in which asymmetric paths are expected to be typical of nonequilibrium steady 
states of interacting particle systems. In Figure 5, the asymmetry seems to grow with N, as in [14|, but a comparison 
of cases with different N is frustrated by the unclear dependence on N of the many parameters entering the definition 
of the asymmetry, like &(X) and frj. For instance, Figure 5 seems to imply that average growth (r < 0) and relaxation 
(t > 0) tend to become linear as N increases, but this may be due to the fact that the observation time to, which ought to 
be connected with the correlations decay rate, increases with N. The right panel of Figure 5 would then only illustrate a 
smaller fraction of the average path than that of the left panel. As a matter of fact, Figure 6 shows that the instantaneous 
asymmetry, 6 T = [p T ~p_ T ], almost saturates with growing t, but it saturates at larger values of r for larger N. On the 
other hand, it is not clear how to should be modified with N; it is only obvious that [p T - p_ T ] must decrease back to 0, 
for t sufficiently long that correlations have decayed JT3J . Unfortunately, that time is too long to allow us to collect a 
statistically relevant sample of fluctuation paths. 



Figure 6. (Color online) Instantaneous asymmetry of average density fluctuations, 6 T = \p T - p_ r ], for different values of the elastic 
constant k 2 , T( = 320, T r = 20 and N = 100 (left panel), N = 200 (central panel), N = 400 (right panel). The asymmetry saturates after 
times which are longer for larger N. 

Let us introduce the asymmetry v, as a normalized cumulative difference between fluctuation and relaxation trajecto- 
ries, 




0.12 




(10) 
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The probability that this quantity be positive, P(v) > 0, is reported in Table 1 for N = 100, where one observes that, from 
this point of view, the k 2 = case does not differ substantially from the k 2 > cases. 



k 2 
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0.8 


0.5 


0.2 


0.1 





P(v > 0) 


0.5660 


0.5716 


0.5486 


0.5710 


0.5648 


0.5462 


0.5823 



Table 1. Probability of positive asymmetry (the fraction of fluctuation paths with v > 0), for A' = 100, T[ = 320, T r = 20, as a function 
of the elastic constant k 2 . 

The above figures and table suggest that nonequilibrium steady states of identical hard spheres on a line do not behave 
like independent particles, although in equilibrium, when the thermostats are removed, they do. Away from equilibrium, 
particles appear to develop space correlations, which reach the bulk, eventually connecting them to the boundaries, even 
in the absence of springs. Therefore, the observed asymmetries reveal that a form of nonlocality (long range correlations) 
[38, 39 40 55 1 is common in nonequilibrium states of both deterministic and stochastic models. 

Another observable of interest is the heat flux. To define the heat flux in the presence of hard core collisions, we adopt 
the method of planes, developed by Todd, Daivis and Evans ||56l . according to which the fluctuating heat flow J q , through 
a plane transversal to the medium, can be decomposed as the sum of two terms: 

J q (x,t) = J*(x,t) + J%(x,i) (11) 

The term j]f represents the kinetic part of the heat flux, which is expressed by: 

1 N 

/£(*, ^ = jYjYj Ui 5{t ~ ^ si Sn[c xi (ti, mi )] (12) 

/= 1 m, 

where i denotes a particle crossing the plane at (discrete) times f, , c A , is the component along direction x of the particle's 
velocity with respect to the streaming velocity, A is the area of the plane (if the system is 3 dimensional), and 

U i = T ^[v i -u{xdf + \<pij (13) 

is the internal energy concerning particle i. Moreover, v, is the velocity of the particle in the laboratory frame, u(xi) is the 
streaming velocity at position x,, and cpij is the interaction potential between particles i and j. The term J q represents the 
potential part of the heat flux, which is defined by: 

1 N 

Jqix, f) = -— ^ [Vj - u(xd] Ftj [sign(* - x t ) - sign(* - Xj)] (14) 

where Fy is the force that particle j exerts on particle i. 

Figure 7 shows the average fluctuations of for the cases with = 100, T( = 320, T r = 20, and for different values 
of k 2 . Analogously to the density average path, the average fluctuation rises faster than it relaxes, and is not appreciably 
affected by the value of k 2 . Results similar to those reported have been obtained for all different parameters choices, even 
with k 2 = 0, when the bulk behaviour should have closely resembled that of non-interacting particles. This is consistent 
with the discussion of the single file behaviour of (53), according to which the order of particles introduces persistent 
correlations, even in the presence of positive Lyapunov exponents. The results appear to be robust and independent of the 
time reversal parity of the observable at hand, since local density and heat flux have opposite parity. 

The results of this section definitely prove that asymmetric fluctuation paths are typical of deterministic, time reversal 
invariant nonequilibrium models. 
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Figure 7. (Color online) Average fluctuations of , for various values of k 2 , with N = 100, T( = 320, T, = 20 and 6 { = $ r = 1, 
The fluctuations have been rescaled, so that they all vanish at time -r, and take the peak value 1 at time 0. 
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4 Harmonic chain in a nonequilibrium steady state 



In this section we consider the model studied by Rieder, Lebowitz and Lieb [ 16 1, with the stochastic boundary thermostats 
replaced by Nose-Hoover thermostats, which make the dynamics time reversal invariant. The oscillators are thus coupled 
by the potential energy V of Eq.tHJ, with (3 — 0. In equilibrium, this system is equivalent to a set of independent particles, 
the normal modes. 

If applicable, the local version of the virial theorem now implies [33|: 

ip*) = ~{(Fi-F M ) qi ), i = 2,...,N-\. (15) 

Clearly, in the absence of thermostats, the local virial theorem does apply to each independent mode of oscillation, but 
in general, this has to be checked. As in the case of stochastic baths considered by |fl6l , we find that Eq.(|T3]> holds and 
that the kinetic temperature profile is characterized by jumps at the boundaries, and by approximately flat profiles in the 
bulk. As in (16'|, harmonic chains do not sustain a kinetic temperature gradient but, differently from the case of [ 16 1, the 
approximate equipartition of energy in the bulk does not correspond to the mean of the left and right temperatures: it stays 
closer to the higher temperature, cf. Fig. 8. 
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Figure 8. (Color online) Temperature profile as a function of the thermostat characteristic time 6. The shorter 0, the stronger the 
coupling between particle and thermostat. Results at short 6 = Of = 6 r are reported in the left half of the left panel, while results at 
large 6 are reported in the right half. After some oscillation, decreasing 6 makes the profile settle around the value 240, which is much 
closer to T t = 320 than to T T = 20. Note the apparent discontinuity of the profile between 6 = 0.5 and 6 = 0.3, which leads to the flat 
bulk profile at (T e + T, )/2. The right panel portrays the profiles with { = 1 and decreasing 6 T ; the system equilibrates with the hot 
thermostat independently of 6 r . In both panels N = 100. 



Some dependence of the bulk behaviour on the characteristics of the thermostats is not surprising: the choice of the 
parameters determines the efficiency of the energy exchange between thermostats and thermostatted system. This, in turn, 
may affect the behaviour of a non-macroscopic system since, by definition, this does not enjoy any local thermodynamic 
equilibrium. In our case, the relevant parameters consist of the target temperatures T[ and T n of the thermostat character- 
istic times 9( and 6 r , and of the elastic constant k 2 . We have studied the behaviour of the profile in the = 8( — 6 r — > 
limit, the limit of strong coupling between thermostats and thermostatted particles, as well as in weak coupling cases. 
Figure 8 (left panel), shows that the kinetic temperature profile settles closer to 7V than to T r , at all values of 8, except 
at 6 = 0.3. The situation does not change if 6 r is shorter than 8(, so that the cold thermostat is more strongly coupled to 
the system than the hot thermostat, as illustrated by the right panel of Fig. 8. For 6 = 0.3, the deterministic time reversal 
invariant Nose-Hoover thermostats behave most similarly to the irreversible stochastic ones of lH6l . as far as the kinetic 
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temperature profile is concerned. Note, in particular, the flatness of the 9 = 0.3 profile, and the temperature overshoot on 
the cold side, predicted by the theory of 0161 . The equivalence with the results for stochastic thermostats of lfl6l is not 
complete, however, because the temperature does not overshoot in the hot side. 
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Figure 9. (Color online) Temperature profile as a function of = 6( = 9 r , for N = 100, TV = 320, T, = 300. For high 9. i.e. 
weak coupling, it settles within [T r , TV], while for low 9, i.e. strong coupling, it reaches even lower than T r . The profile depends quite 
irregularly on 9. For 9 = 0.7 and 9=1 the bulk temperature settles close to {T e + T r )/2. 

This is consistent with the observations of Ref.[47|, on a system of rotating disks and pointlike particles, thermostatted 
by Nose-Hoover mechanisms acting on the boundary disks. There, the hotter side reaches local equilibrium more easily 
than the colder one, even if not necessarily at the target temperature of the thermostat. The thermostat closer to an 
equilibrium state (as it should be) seems, then, more efficient in driving the bulk of the chain than the other thermostat. In 
Ref. BTl , it was further observed that small gradients lead to standard behaviour of the cold thermostat. Differently, in our 
case, the temperature profiles remain quite irregular even under relatively small gradients, and only for a limited set of 9's 
do they lie close to (7> + T r )/2. The profiles worsen in the small 9 limit, i.e. for strong couplings, when they even brim 
over the interval [T r , 7V], cf. Fig. 9. The strong dependence on the microscopic parameters is confirmed by the departure 
at large 9 of the bulk profiles from the mean of the boundary values, cf. the case 9 — 2 in the left panel of Fig. 9. 

To further assess the relation between our chains and those of [ 16] we have produced the histograms of pi, and have 
verified that they are not Gaussian, although to different degrees for the different i's. Indeed, the distributions of pi 
and pn/2 differ much less from a normal distribution than that of p^ (N is the particle interacting with the cold bath 
at temperature T[), as indicated by normal distribution plots not reported here. Also, local Gaussian distributions are 
better and better approximated at both ends of the chain, as in [47 1, if the temperature gradient is reduced. Increasing k 2 , 
without changing the thermostats properties, is another route towards approximately Gaussians local distributions. This 
confirms that closer to equilibrium it is easier to obtain the equivalence between deterministic and stochastic thermostats, 
as generally expected (26], although the temperature profiles demonstrate that complete equivalence is not to be expected. 

The peculiar behaviour of the purely harmonic nonequilibrium chains is evidenced also by quantities such as (xi+\ - xi) 
and {qi), whose profiles are much more irregular than those of Figure 3, cf. Figure 10. Here, differently from Figure 3, 
(qt) almost vanishes, consistently with the approximately equilibrium state of the bulk of the chain@ 

'Note that this state only looks close to equilibrium: for instance, position- velocity correlations differ significantly from their equilibrium counterparts 
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ure 10. (Color online) Average coordinates difference between nearest neighbours, {xi+\ - x,-), for k 2 = 1 and equilibrium distance 
5 (left panel). Average displacement (qt) from the equilibrium position, ai (right panel). In both panels N = 100, T( = 320, 
= 20, k 1 = 1. 
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5 Lyapunov exponents: the harmonic and the hard core cases 



In this section we consider the Lyapunov exponents of the two limiting cases of purely harmonic and purely hard core in- 
teractions. These exponents have been computed in tangent space, with the usual algorithm devised by Benettin, Galgani, 
Giorgilli and Strelcyn 11571 . The phase space contraction rate^, i.e. the dynamical dissipation, can be separately computed 
as the negative of the average of the divergence of the vector field, neglecting the instantaneous elastic collisions, which 
do not contribute to the variations of the phase space volumes |58|. In our case, (x) = ((&) + {&)) and (x) > 0, where 
equality characterizes equilibrium states, while (y) > characterizes nonequilibrium states. 

The motion of purely harmonic chains, in the absence of thermostats, is hamiltonian and fully integrable, hence all 
Lyapunov exponents vanish. Table 2 shows the values of the largest Lyapunov exponent, A\, for several different kinetic 
temperatures at the left end of a chain of 100 particles, while the right end is subjected to a thermostat with T r - 20. 



Tt 


M 


(x) 


420 


0.00771 ±0.5 x 10~ 4 


2.7988 


370 


0.00769 +0.6 x 10~ 4 


2.5737 


320 


0.00779 +0.7 x 10~ 4 


2.3382 


220 


0.00792 +0.5 x 10~ 4 


1.7644 


120 


0.00928+0.1 x 10~ j 


1.0726 


80 


0.01101 +0.1 x 10" a 


0.7194 


50 


0.01338 +0.2 x 10~ j 


0.3775 


20 


0.01328 +0.1 x 10~ J 


6.54 x 10 -5 



Table 2. Largest Lyapunov exponent, Ai, and average phase space contraction rate, (x), as functions of T ( , for purely harmonic chains 
with A' = 100, k 1 = 1 and T, = 20. Standard deviations of the computed values are reported, with the number in brackets representing 
the power of 10 which multiplies them. 

It is interesting to note that even at equilibrium, i.e. for Tg = T r , the presence of the thermostats makes chaotic the 
dynamics. Indeed, A\ is higher at lower temperature differences, since higher dissipations imply smaller attractors and 
more orderly states. It also appears that A\ saturates as a function of the temperature difference, while the dissipation keeps 
increasing significantly. As a function of the relaxation times of the thermostats, the largest Lyapunov exponent decreases 
mono tonic ally. For instance, the largest exponent for N = 100 and T{ = 320 decreases from the value A\ = 0.00779 
reported in Table 2 for = 1 to A\ = 0.000444 for 9 = 10. This is a consequence of the fact that the motion of single 
oscillators is integrable and that large 9 practically decouples the thermostats from the system of interest. Table 3 reports 
the computed values of A\ and x, at fixed temperature difference and varying chain length. The Lyapunov exponent 
decreases with N, hence with the temperature gradient, while the dissipation (x) is practically constant. 



N 


Ay 


<x) 


50 


0.0120 ±0.8 x 10~ 4 


2.3013 


100 


0.00779 ±0.7 x 10~ 4 


2.3382 


150 


0.0062 ±0.8 x 10~ 4 


2.3286 


200 


0.0048 ±0.1 x 10~ J 


2.3296 



Table 3. Largest Lyapunov exponent and dissipation for the purely harmonic case, with k 2 = 1, and fixed temperature difference: 
T r = 20, T[ = 320. The global gradient decreases with N, like At, while (x) is practically constant. 

The picture is completed by the behaviour with N at fixed temperature gradient, cf. Table 40 The data indicate that 
two competing effects contribute to the values of the largest Lyapunov exponent: 

10 As observed in e.g. [59], the local kinetic temperature gradient inside the chain is not a useful variable, in this case, because it fluctuates. 
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a) decreasing the temperature gradient, the system approaches the equilibrium state, which enjoys the largest value of 
Ai afforded by a chain of given length; 

b) increasing N reduces the impact of the thermostats on the bulk of the chain, which then better approximates an 
isolated, fully integrable system, with vanishing Lyapunov exponents. 

As Ai does not grow indefinitely with (7> - T r ), increasing N leads to a decrease of A\. The dissipation {%), on the other 
hand, is more directly related to the value of (Tf - 7V), and grows with it. 



N 


M 


(X) 


50 


0.0131 +0.68 x 10- 4 


1.4384 


100 


0.0078 ±0.73 x 10~ 4 


2.3382 


150 


0.0057 +0.96 x 10~ 4 


3.0205 


200 


0.0050 +0.67 x 10~ 4 


3.6501 



Table 4. Largest Lyapunov exponent and dissipation for the purely harmonic case, with k 1 = 1, T r = 20 and fixed global temperature 
gradient (7> - T r )/N = 3. 

An interesting feature of the purely harmonic case, with Nose-Hoover thermostats at different temperatures, is that 
non-chaotic dynamics can at least be realized in short chains. For instance, in the case with N = 3, T t = 80, k 2 = 1 and 
T r = 20, all Lyapunov exponents are negative, except for one vanishing exponent, and the motion is periodic. As a matter 
of fact, even some cases with larger N might share these features, but their largest exponents appear to be only marginally 
positive, hence hard to distinguish from vanishing values. Differently, the case with N - A has clearly positive exponents. 

Consider now the case in which particles interact via hard core elastic collisions only, which, in the absence of 
thermostats, would be the other limiting, fully integrable case. The largest exponents are commonly found to be positive 
and behave as in the purely harmonic case, cf. Table 5. 



T t 


M 


<x) 


320 


0.01641+0.16 x \Q- 2 


8.41388 


220 


0.01931+0.18 x 10^ 


4.93916 


120 


0.01978+0.10 x 10~ J 


2.05788 


80 


0.02441 ±0.14 x 10~ J 


1.02030 


50 


0.02983 ±0.22 x 10~ J 


0.33846 


20 


0.03157 ±0.16 x 10~ J 


0.04161 



Table 5. Largest Lyapunov exponent and dynamical dissipation for hard particles with k 2 = 0, N = 100 and T r = 20. The behaviours 
of both A[ and (%) do not differ substantially from those of purely harmonic cases. 

The comparison between hard core and purely harmonic cases illustrates one interesting fact, consistent with the ob- 
servations of Ref. 1411 . Both kinds of systems have positive Lyapunov exponents, but the purely harmonic case shows a 
quite peculiar bulk equipartition of energy, while the chains of purely hard core particles verify more standard nonequi- 
librium conditions. The randomness entailed by hard collisions, which occur at practically random times and positions, 
as assumed in [27], breaks correlations among particles more efficiently than generic chaotic mechanisms do and favours 
standard behaviour, as argued also in [41 1. However, the other peculiarities noted in the previous sections still do not 
allow us to speak of genuine local thermodynamic equilibrium for chains of hard particles. That chaos does not suffice 
for standard behaviour, in FPU systems, has been observed in the past 13311431 . 
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6 Concluding remarks 



We have analyzed deterministic chains of oscillators with harmonic forces and elastic collisions, to asses the equivalence 
of deterministic time reversal invariant and stochastic models of thermostats. Close to equilibrium, it is easier to obtain 
equivalent behaviours than far from equilibrium, although the peculiarities of 1 -d dynamics prevent complete equivalence, 
even in the large N limit. Only from a qualitative standpoint, various properties are common to the different models, like 
the form of the temperature profiles, the temporal asymmetries of flucutations and the nonlocality of nonequilibrium steady 
states [3,8, 9 . 38 40 1. Similarly, the strong dependence of the bulk behaviour on microscopic details of the dynamics, such 
as the boundary conditions, which is enhanced by the growth of the dissipation, reveals qualitative similarities between 
the different kinds of thermostatted evolutions. This qualitative similarity is robust against changes of the interaction 
potentials. For instance, the overall behaviour of systems made of purely harmonic identical oscillators differs from that 
implied by hard core and other non-linear forces, but this happens for both stochastic and deterministic thermostats. 

Therefore, from relatively coarse, qualitative, viewpoints may we state that the different models of thermostat are 
equivalent. This weak equivalence suffices to conclude that anomalous behaviours are in fact typical of physics in (quasi) 
1 -dimension, but does not suffice to predict the kind of anomalies that one should expect in different situations. As a 
matter of fact, Appendix B shows that numerous parameters affect the temperature profiles, including the elastic constant 
k 2 and any perturbation of the interaction forces. In particular, the presence of different degrees of stochasticity produces 
different behaviours, and only in the equilibrium limit may the equivalence of the different models be established. 

Independently of the nature of the thermostats, "realistic" particle interactions, or mechanisms preventing the conser- 
vation of the most obvious dynamical quantities do not lead to the onset of local thermodynamic equilibrium l|26l I4T1 |46l . 
One reason is that certain correlations (e.g. particles order) persist in time, violating the molecular chaos hypothesis of 
kinetic theory, in 1-d or quasi- 1-d systems. Thus, some kind of anomalous behaviour is to be generically expected, as a 
consequence of the low dimensionality of the dynamics, rather than of the peculiarities of thermostats. Indeed, in 3 di- 
mensions, the same thermostatting mechanisms do not lead to anomalous behaviour. Even in cases of apparently normal 
transport in 1-d, this is not as robust, with respect to variations of the microscopic parameters, as it is in macroscopic ther- 
modynamic phenomena ll26ll4TI . indicating that genuine local thermodynamic equilibrium does not hold in those cases. 
Experimental evidence, although presently scanty, confirms this picture, cf. e.g. 1461151 52 ,53 ,60, 61 , 62 1. Further study 
is desired, to understand the relation between the peculiarities of the 1-d models so far considered in the literature and 
those of real systems, not in local thermodynamic equilibrium. 

We found the linear relation (0 between the particles mean displacements from their equilibrium positions and the 
temperature profiles. Apart from the temporal asymmetries of the fluctuations of the main observables, considered in 
Section 3, relation (0 is the only result which does not show a delicate dependence on the details of the microscopic 
dynamics. It is robust against all modifications of the dynamics, which we have investigated, and must therefore play an 
important role in the identification of universality classes and in the equivalence of thermostats in 1-d systems 

Our analysis of systems with purely harmonic and purely hard core interactions indicates that chaos, or a generic 
source of randomness, per se, does not suffice to establish standard nonequilibrium steady states. On the other hand, it 
has been observed by various authors that dynamical chaos is not even necessary for that, see e.g. ll4Tll42l . 
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Appendix A 



To assess the relevance of our results for the large N limit, and to compare with the results of Refs. 11271 l28l we have 
produced temperature profiles for hard core collisions and k 2 = 1, with TV = 320, T r = 20, Qt — 8 r — 1, and growing N. 
Figure 1 1 shows these profiles, plotted as functions of the rescaled variable x = i/N. One observes a scaling similar to 
that of Fig. 1 in Ref.[43], where the asymptotic profile is well approximated, on the scale of the figure, at moderately large 
values of N. 




Magnifying this figure, one observes that the profile slowest convergence occurs close to the cold boundary, cf. Fig. 12. 
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Figure 12. (Color online) Magnification of Fig. 1 1 concerning the center of the chain (left panel) and the cold side (right panel). 



The difference AjvCy) = T(y) - T^(y) between the analytic profile T of Ref.[27| and our profiles, 7V, with normalized 
variable ye [-1,1], for finite N, does not scale with as the difference between asymptotic and finite profiles of ll27l . 
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If the scaling was confirmed, we would have had 



AatOO oc N- l/i g(y) (16) 

where g(y) depends on the details of the model, but not on N. As a matter of fact, the absolute value of A N (y) initially 
grows with N, rather than decreasing, e.g. in the center of chain. This does not prevent the scaling to set in at still larger 
values of N, although this looks unlikely, given the observed convergence of our profiles to a different asymptotic shape. 



20 



Appendix B 



We consider a stochastic perturbation of the model of Section 2. If the dynamics of Section 2 leads two particles to collide 
at a given instant of time f, they will actually collide with probability 1 - p, p e [0, 1], at time t. Taking p = yields 
the dynamics of Section 2, while p — 1 yields the purely harmonic dynamics of Section 4, because no collision takes 
place. For a system of N - 100 particles, we find that the temperature profiles strongly depend on the values of the pair 
of parameters p and k 2 , confirming that nonequilibrium chains of oscillators may hardly be part of a unique universality 
class, except in the equilibrium limit, i.e. for temperatures at the boundaries very close to each other. The states of these 
chains are indeed too far from local equilibrium, and sensitive to many details of the microscopic dynamics, in general. 

In particular, Fig. 13 shows the temperature profiles for different values of p and k 2 = 0.1, for two temperature 
differences (T(-T r ). The dependence on p is quite strong although, in both cases, the profiles settle around the unperturbed 
ones, when p decreases. The cases with lower temperature difference more rapidly collapse on a unique shape, different 
from the theoretical T of [1271 . 





Figure 13. (Color online) Temperature profiles as functions of p, for N = 100, k 2 =0.1. The left panel reports the case with T( = 320 
and T r = 20. The right panel reports the case with 77 = 120 and T, = 20. The continuous line represents the temperature profile T(x) 
of eq. ([SJ. The lower temperature jump leads to more rapid convergence to the unperturbed profiles, as p — > 0. 



Figure 14, shows that increasing the rigidity of the chains, reduces the effect of the parameter p. Comparing Figs. 13 
and 14 leads to the conclusion that, even for a fixed kind of thermostats, the equivalence of the different microscopic 
dynamics requires various parameters to be adjusted, in general. In our case, k 2 plays an important role, which it did not 
in Refs.[27 . 28 1. Taking of k 2 — 2 does not yield profiles which approximate the theoretical one of [27] more closely than 
^=1. 



21 




Figure 14. (Color online) Temperature profiles as functions of p, for N = 100, k 2 = 1. The left panel reports the case with T( = 320 
and T, = 20. The right panel reports the case with TV = 120 and T r = 20. The continuous line represents the temperature profile T(x) 
of eq. l|6}. The dependence on p is reduced at this higher value of k 2 . 



The only result which is robust against all the modifications of the dynamics, which we have investigated, is the validity 
of the linear relation (0, which must then play an important role in the identification of universality classes and in 
the equivalence of thermostats in 1-d systems. Figures 15 and 16 illustrate this fact, confirmed by all simulations we 
performed, for only two of these cases. 




Figure 15. (Color online) Mean displacement (<7;)t2 = o.i from the equilibrium centres of oscillation of the particles, ai (left panel) and 
mean displacement {qi) k i =l (right panel), for N = 100, T ( = 320 and T T = 20. Larger k 2 reduces the p dependence, except for p = 0, 
which does not sustain temperature gradients. 



For a given choice of the parameters of the deterministic model, the values of the parameters j3\ and of Eq.© do not 
depend on p, as long as p < 0.9. They change discontinuosly for p > 0.9, i.e. close to the purely harmonic chains, which 
are quite peculiar models, compared to the others. 
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50 i 100 50 i 100 

Figure 16. (Color online) Relation between mean coordinates difference i - x,) and temperature profiles, for N = 100, 7> = 320, 
T, = 20, k 2 = 0.1. The left panel reports the case p = 0.1, the right panel reports the case p = 0.9. Relation l[5} is confirmed in all 
instances as discussed in Section 2. 
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